Conformational Tuning Shapes the Balance between Functional Promiscuity and Specialization in Paralogous Plasmodium Acyl-CoA Binding Proteins.

Paralogous proteins confer enhanced fitness to organisms via complex sequence-conformation codes that shape functional divergence, specialization, or promiscuity. Here, we dissect the underlying mechanism of promiscuous binding versus partial subfunctionalization in paralogues by studying structurally identical acyl-CoA binding proteins (ACBPs) from Plasmodium falciparum that serve as promising drug targets due to their high expression during the protozoan proliferative phase. Combining spectroscopic measurements, solution NMR, SPR, and simulations on two of the paralogues, A16 and A749, we show that minor sequence differences shape nearly every local and global conformational feature. A749 displays a broader and heterogeneous native ensemble, weaker thermodynamic coupling and cooperativity, enhanced fluctuations, and a larger binding pocket volume compared to A16. Site-specific tryptophan probes signal a graded reduction in the sampling of substates in the holo form, which is particularly apparent in A749. The paralogues exhibit a spectrum of binding affinities to different acyl-CoAs with A749, the more promiscuous and hence the likely ancestor, binding 1000-fold stronger to lauroyl-CoA under physiological conditions. We thus demonstrate how minor sequence changes modulate the extent of long-range interactions and dynamics, effectively contributing to the molecular evolution of contrasting functional repertoires in paralogues.

Roles of acyl-CoA-binding proteins in plant reproduction.

Acyl-CoA-binding proteins (ACBPs) constitute a well-conserved family of proteins in eukaryotes that are important in stress responses and development. Past studies have shown that ACBPs are involved in maintaining, transporting and protecting acyl-CoA esters during lipid biosynthesis in plants, mammals, and yeast. ACBPs show differential expression and various binding affinities for acyl-CoA esters. Hence, ACBPs can play a crucial part in maintaining lipid homeostasis. This review summarizes the functions of ACBPs during the stages of reproduction in plants and other organisms. A comprehensive understanding on the roles of ACBPs during plant reproduction may lead to opportunities in crop improvement in agriculture.

Investigations of Lipid Binding to Acyl-CoA-Binding Proteins (ACBP) Using Isothermal Titration Calorimetry (ITC).

Isothermal titration calorimetry (ITC) is a quantitative, biophysical method to investigate intermolecular binding between biomolecules by directly measuring the heat exchange in the binding reaction. The assay is carried out in solution when the molecules interact in vitro. This allows to determine values for binding affinity (Kd), binding stoichiometry (n), as well as changes in Gibbs free energy (ΔG), entropy (ΔS), and enthalpy (ΔH). This method also addresses the kinetics of enzymatic reactions for a substrate during conversion to a product. ITC has been used to study the interactions between proteins and ligands such as those of acyl-CoA-binding proteins (ACBPs) and acyl-CoA thioesters or ACBPs with protein partners. ITC has also been used in investigating interactions between antiserum and antigen, as well as those involving RNA and DNA and other macromolecules. We describe the methods used to isolate and purify a recombinant rice ACBP (OsACBP) for ITC. To study OsACBP binding to long-chain acyl-CoA thioesters, a microcalorimeter was used at 30 °C, and the ligand (acyl-CoA thioesters or a protein partner in the first cell), was mixed with the ACBP protein solution in a second cell, for more than 40 min comprising 20 injections. Subsequently, the binding parameters including the heat-release data were analyzed and various thermodynamic parameters were calculated.

Robust and Accurate Computational Estimation of the Polarizability Tensors of Macromolecules.

Alignment of molecules through electric fields minimizes the averaging over orientations, e.g., in single-particle-imaging experiments. The response of molecules to external ac electric fields is governed by their polarizability tensor, which is usually calculated using quantum chemistry methods. These methods are not feasible for large molecules. Here, we calculate the polarizability tensor of proteins using a regression model that correlates the polarizabilities of the 20 amino acids with perfect conductors of the same shape. The dielectric constant of the molecules could be estimated from the slope of the regression line based on the Clausius-Mossotti equation. We benchmark our predictions against the quantum chemistry results for the Trp cagemini protein and the measured dielectric constants of larger proteins. Our method has applications in computing laser alignment of macromolecules, for instance, benefiting single-particle imaging, as well as for estimation of the optical and electrostatic characteristics of proteins and other macromolecules.

Exploring the Denatured State Ensemble by Single-Molecule Chemo-Mechanical Unfolding: The Effect of Force, Temperature, and Urea.

While it is widely appreciated that the denatured state of a protein is a heterogeneous conformational ensemble, there is still debate over how this ensemble changes with environmental conditions. Here, we use single-molecule chemo-mechanical unfolding, which combines force and urea using the optical tweezers, together with traditional protein unfolding studies to explore how perturbants commonly used to unfold proteins (urea, force, and temperature) affect the denatured-state ensemble. We compare the urea m-values, which report on the change in solvent accessible surface area for unfolding, to probe the denatured state as a function of force, temperature, and urea. We find that while the urea- and force-induced denatured states expose similar amounts of surface area, the denatured state at high temperature and low urea concentration is more compact. To disentangle these two effects, we use destabilizing mutations that shift the Tm and Cm. We find that the compaction of the denatured state is related to changing temperature as the different variants of acyl-coenzyme A binding protein have similar m-values when they are at the same temperature but different urea concentration. These results have important implications for protein folding and stability under different environmental conditions.

Analysis of an acyl-CoA binding protein in Aspergillus oryzae that undergoes unconventional secretion.

Acyl-CoA binding protein (ACBP) plays important roles in the metabolism of lipids in eukaryotic cells. In the industrially important filamentous fungus Aspergillus oryzae, although we have previously demonstrated that the A. oryzae ACBP (AoACBP) localizes to punctate structures and exhibits long-range motility, which is dependent on autophagy-related proteins, the physiological role of AoACBP remains elusive. Here, we describe identification and characterization of another ACBP from A. oryzae; we named this ACBP as AoAcb2 and accordingly renamed AoACBP as AoAcb1. The deduced amino acid sequence of AoAcb2 lacked a signal peptide. Phylogenetic analysis classified AoAcb2 into a clade that was same as the ACBP Acb1 of the model yeast Saccharomyces cerevisiae, but was different from that of AoAcb1. In contrast to punctate localization of AoAcb1, AoAcb2 was found to be dispersedly distributed in the cytoplasm, as was previously observed for the S. cerevisiae Acb1. Since we could not generate an Aoacb2 disruptant, we created an Aoacb2 conditional mutant that exhibited less growth under Aoacb2-repressed condition, suggesting that Aoacb2 is an essential gene for growth. Moreover, we observed that A. oryzae AoAcb2, but not A. oryzae AoAcb1, was secreted under carbon-starved condition, suggesting that AoAcb2 might be secreted via the unconventional protein secretion (UPS) pathway, just like S. cerevisiae Acb1. We also demonstrated that the unconventional secretion of AoAcb2 was dependent on the t-SNARE AoSso1, but was independent of the autophagy-related protein AoAtg1, suggesting that the unconventional secretion of AoAcb2, unlike that of S. cerevisiae Acb1, via the UPS pathway, is not regulated by the autophagy machinery. Thus, the filamentous fungus A. oryzae harbors two types of ACBPs, one of which appears to be essential for growth and undergoes unconventional secretion.

Lipid membranes and acyl-CoA esters promote opposing effects on acyl-CoA binding protein structure and stability.

Acyl-CoA Binding Proteins (ACBP) form a housekeeping family of proteins that is responsible for the buffering of long chain acyl-coenzyme A esters (LCFA-CoA) inside the cell. Even though numerous studies have focused on the characterization of different members of the ACBP family, the knowledge about the impact of both LCFA-CoA and phospholipids on ACBP structure and stability remains scarce. Besides, there are still controversies regarding the possible interaction of ACBP with biological membranes, even though this might be essential for the cargo capture and delivery. In this study, we observed that LCFA-CoA and phospholipids play opposite roles on protein stability and that the interaction with the membrane is dictated by electrostatic interaction. Furthermore, the results support the hypothesis that the LCFA-CoA delivery is driven by the increase of the negative charge on the membrane surface. The combined influence played by the different molecules on ACBP structure is discussed on the light of cargo capture/delivery giving new insights about this important process.

Corneal haze phenotype in Aldh3a1-null mice: In vivo confocal microscopy and tissue imaging mass spectrometry.

ALDH3A1 is a corneal crystallin that protects ocular tissues from ultraviolet radiation through catalytic and non-catalytic functions. In addition, ALDH3A1 plays a functional role in corneal epithelial homeostasis by simultaneously modulating proliferation and differentiation. We have previously shown that Aldh3a1 knockout mice in a C57B6/129sV mixed genetic background develop lens cataracts. In the current study, we evaluated the corneal phenotype of Aldh3a1 knockout mice bred into a C57B/6J congenic background (KO). In vivo confocal microscopy examination of KO and wild-type (WT) corneas revealed KO mice to exhibit corneal haze, manifesting marked light scattering from corneal stroma. This corneal phenotype was further characterized by Imaging Mass Spectrometry (IMS) with spatial resolution that revealed a trilayer structure based on differential lipid localization. In these preliminary studies, no differences were observed in lipid profiles from KO relative to WT mice; however, changes in protein profiles of acyl-CoA binding protein (m/z 9966) and histone H4.4 (m/z 11308) were found to be increased in the corneal epithelial layer of KO mice. This is the first study to use IMS to characterize endogenous proteins and lipids in corneal tissue and to molecularly explore the corneal haze phenotype. Taken together, the current study presents the first genetic animal model of cellular-induced corneal haze due to the loss of a corneal crystallin, and strongly supports the notion that ALDH3A1 is critical to cellular transparency. Finally, IMS represents a valuable new approach to reveal molecular changes underlying corneal disease.

Recombinant expression, purification, and characterization of an acyl-CoA binding protein from Aspergillus oryzae.

OBJECTIVES: To characterize biochemically the lipid metabolism-regulating acyl-CoA binding protein (ACBP) from the industrially-important fungus Aspergillus oryzae. RESULTS: A full-length cDNA encoding a candidate ACBP from A. oryzae (AoACBP) was cloned and expressed in Escherichia coli as a maltose-binding protein (MBP) fusion protein. The MBP-AoACBP protein was purified by an amylose resin chromatography column. SDS-PAGE showed that MBP-AoACBP has an estimated molecular weight of 82 kDa. Microscale thermophoresis binding assay showed that the recombinant AoACBP displayed much greater affinity for palmitoyl-CoA (K d = 80 nM) than for myristoyl-CoA (K d = 510 nM), thus demonstrating the preference of AoACBP for long-chain acyl-CoA. CONCLUSION: The data support the identification of AoACBP as a long-chain ACBP in A. oryzae.

Genome-wide identification and Phylogenic analysis of kelch motif containing ACBP in Brassica napus.

BACKGROUND: Acyl-coA binding proteins (ACBPs) bind long chain acyl-CoA esters with very high affinity. Their possible involvement in fatty acid transportation from the plastid to the endoplasmic reticulum, prior to the formation of triacylglycerol has been suggested. Four classes of ACBPs were identified in Arabidopsis thaliana: the small ACBPs, the large ACBPs, the ankyrin repeats containing ACBPs and the kelch motif containing ACBPs. They differed in structure and in size, and showed multiple important functions. In the present study, Brassica napus ACBPs were identified and characterized. RESULTS: Eight copies of kelch motif ACBPs were cloned, it showed that B. napus ACBPs shared high amino acid sequence identity with A. thaliana, Brassica rapa and Brassica oleracea. Furthermore, phylogeny based on domain structure and comparison map showed the relationship and the evolution of ACBPs within Brassicaceae family: ACBPs evolved into four separate classes with different structure. Chromosome locations comparison showed conserved syntenic blocks. CONCLUSIONS: ACBPs were highly conserved in Brassicaceae. They evolved from a common ancestor, but domain duplication and rearrangement might separate them into four distinct classes, with different structure and functions. Otherwise, B. napus inherited kelch motif ACBPs from ancestor conserving chromosomal location, emphasizing preserved synteny block region. This study provided a first insight for exploring ACBPs in B. napus, which supplies a valuable tool for crop improvement in agriculture.

Computational Prediction of acyl-coA Binding Proteins Structure in Brassica napus.

Acyl-coA binding proteins could transport acyl-coA esters from plastid to endoplasmic reticulum, prior to fatty acid biosynthesis, leading to the formation of triacylglycerol. The structure and the subcellular localization of acyl-coA binding proteins (ACBP) in Brassica napus were computationally predicted in this study. Earlier, the structure analysis of ACBPs was limited to the small ACBPs, the current study focused on all four classes of ACBPs. Physicochemical parameters including the size and the length, the intron-exon structure, the isoelectric point, the hydrophobicity, and the amino acid composition were studied. Furthermore, identification of conserved residues and conserved domains were carried out. Secondary structure and tertiary structure of ACBPs were also studied. Finally, subcellular localization of ACBPs was predicted. The findings indicated that the physicochemical parameters and subcellular localizations of ACBPs in Brassica napus were identical to Arabidopsis thaliana. Conserved domain analysis indicated that ACBPs contain two or three kelch domains that belong to different families. Identical residues in acyl-coA binding domains corresponded to eight amino acid residues in all ACBPs of B. napus. However, conserved residues of common ACBPs in all species of animal, plant, bacteria and fungi were only inclusive in small ACBPs. Alpha-helixes were displayed and conserved in all the acyl-coA binding domains, representing almost the half of the protein structure. The findings confirm high similarities in ACBPs between A. thaliana and B. napus, they might share the same functions but loss or gain might be possible.

Imprints of function on the folding landscape: functional role for an intermediate in a conserved eukaryotic binding protein.

In the computational characterization of single domain protein folding, the effective free energies of numerous microstates are projected onto few collective degrees of freedom that in turn serve as well-defined reaction coordinates. In this regard, one-dimensional (1D) free energy profiles are widely used mainly for their simplicity. Since folding and functional landscapes are interlinked, how well can these reduced representations capture the structural and dynamic features of functional states while being simultaneously consistent with experimental observables? We investigate this issue by characterizing the folding of the four-helix bundle bovine acyl-CoA binding protein (bACBP), which exhibits complex equilibrium and kinetic behaviours, employing an Ising-like statistical mechanical model and molecular simulations. We show that the features of the 1D free energy profile are sufficient to quantitatively reproduce multiple experimental observations including millisecond chevron-like kinetics and temperature dependence, a microsecond fast phase, barrier heights, unfolded state movements, the intermediate structure and average ϕ-values. Importantly, we find that the structural features of the native-like intermediate (partial disorder in helix 1) are intricately linked to a unique interplay between packing and electrostatics in this domain. By comparison with available experimental data, we propose that this intermediate determines the promiscuous functional behaviour of bACBP that exhibits broad substrate specificity. Our results present evidence to the possibility of employing the statistical mechanical model and the resulting 1D free energy profile to not just understand folding mechanisms but to even extract features of functionally relevant states and their energetic origins.

Exploring the minimally frustrated energy landscape of unfolded ACBP.

The unfolded state of globular proteins is not well described by a simple statistical coil due to residual structural features, such as secondary structure or transiently formed long-range contacts. The principle of minimal frustration predicts that the unfolded ensemble is biased toward productive regions in the conformational space determined by the native structure. Transient long-range contacts, both native-like and non-native-like, have previously been shown to be present in the unfolded state of the four-helix-bundle protein acyl co-enzyme binding protein (ACBP) as seen from both perturbations in nuclear magnetic resonance (NMR) chemical shifts and structural ensembles generated from NMR paramagnetic relaxation data. To study the nature of the contacts in detail, we used paramagnetic NMR relaxation enhancements, in combination with single-point mutations, to obtain distance constraints for the acid-unfolded ensemble of ACBP. We show that, even in the acid-unfolded state, long-range contacts are specific in nature and single-point mutations affect the free-energy landscape of the unfolded protein. Using this approach, we were able to map out concerted, interconnected, and productive long-range contacts. The correlation between the native-state stability and compactness of the denatured state provides further evidence for native-like contact formation in the denatured state. Overall, these results imply that, even in the earliest stages of folding, ACBP dynamics are governed by native-like contacts on a minimally frustrated energy landscape.

Molecular properties of the class III subfamily of acyl-coenyzme A binding proteins from tung tree (Vernicia fordii).

Acyl-CoA binding proteins (ACBPs) have been identified in most branches of life, and play various roles in lipid metabolism, among other functions. Plants contain multiple classes of ACBP genes. The most diverse group is the class III proteins. Tung tree (Vernicia fordii) contains two such genes, designated VfACBP3A and VfACBP3B. The two proteins are significantly different in length and sequence. Analysis of tung ACBP3 genes revealed significant evolution, suggesting relatively ancient divergence of the two genes from a common ancestor. Phylogenetic comparisons of multiple plant class III proteins suggest that this group is the most evolutionarily dynamic class of ACBP. Both tung ACBP3 genes are expressed at similar levels in most tissues tested, but ACBP3A is stronger in leaves. Three-dimensional modeling predictions confirmed the presence of the conserved four α-helix bundle acyl-CoA binding (ACB); however, other regions of these proteins likely fold much differently. Acyl-CoA binding assays revealed different affinities for different acyl-CoAs, possibly contradicting the redundancy of function suggested by the gene expression studies. Subcellular targeting of transiently-expressed plant ACBP3 proteins contradicted earlier studies, and suggested that at least some class III ACBPs may be predominantly targeted to endoplasmic reticulum membranes, with little or no targeting to the apoplast.

Ratcheted molecular-dynamics simulations identify efficiently the transition state of protein folding.

The atomistic characterization of the transition state (TS) is a fundamental step to improve the understanding of the folding mechanism and the function of proteins. From a computational point of view, the identification of the conformations that build out the transition state is particularly cumbersome, mainly because of the large computational cost of generating a statistically sound set of folding trajectories. Here we show that a biasing algorithm, based on the physics of the ratchet-and-pawl, can be used to approximate efficiently the transition state. The basic idea is that the algorithmic ratchet exerts a force on the protein when it is climbing the free-energy barrier, while it is inactive when it is descending. The transition state can be identified as the point of the trajectory where the ratchet changes regime. Besides discussing this strategy in general terms, we test it within a protein model whose transition state can be studied independently by plain molecular dynamics simulations. Finally, we show its power in explicit-solvent simulations, obtaining and characterizing a set of transition-state conformations for Acyl-Coenzyme A-Binding Protein (ACBP) and Chymotrypsin Inhibitor 2 (CI2).

A highly compliant protein native state with a spontaneous-like mechanical unfolding pathway.

The mechanical properties of proteins and their force-induced structural changes play key roles in many biological processes. Previous studies have shown that natively folded proteins are brittle under tension, unfolding after small mechanical deformations, while partially folded intermediate states, such as molten globules, are compliant and can deform elastically a great amount before crossing the transition state barrier. Moreover, under tension proteins appear to unfold through a different sequence of events than during spontaneous unfolding. Here, we describe the response to force of the four-α-helix acyl-CoA binding protein (ACBP) in the low-force regime using optical tweezers and ratcheted molecular dynamics simulations. The results of our studies reveal an unprecedented mechanical behavior of a natively folded protein. ACBP displays an atypical compliance along two nearly orthogonal pulling axes, with transition states located almost halfway between the unfolded and folded states. Surprisingly, the deformability of ACBP is greater than that observed for the highly pliant molten globule intermediate states. Furthermore, when manipulated from the N- and C-termini, ACBP unfolds by populating a transition state that resembles that observed during chemical denaturation, both for structure and position along the reaction coordinate. Our data provide the first experimental evidence of a spontaneous-like mechanical unfolding pathway of a protein. The mechanical behavior of ACBP is discussed in terms of topology and helix propensity.

Slow unfolded-state structuring in Acyl-CoA binding protein folding revealed by simulation and experiment.

Protein folding is a fundamental process in biology, key to understanding many human diseases. Experimentally, proteins often appear to fold via simple two- or three-state mechanisms involving mainly native-state interactions, yet recent network models built from atomistic simulations of small proteins suggest the existence of many possible metastable states and folding pathways. We reconcile these two pictures in a combined experimental and simulation study of acyl-coenzyme A binding protein (ACBP), a two-state folder (folding time ~10 ms) exhibiting residual unfolded-state structure, and a putative early folding intermediate. Using single-molecule FRET in conjunction with side-chain mutagenesis, we first demonstrate that the denatured state of ACBP at near-zero denaturant is unusually compact and enriched in long-range structure that can be perturbed by discrete hydrophobic core mutations. We then employ ultrafast laminar-flow mixing experiments to study the folding kinetics of ACBP on the microsecond time scale. These studies, along with Trp-Cys quenching measurements of unfolded-state dynamics, suggest that unfolded-state structure forms on a surprisingly slow (~100 µs) time scale, and that sequence mutations strikingly perturb both time-resolved and equilibrium smFRET measurements in a similar way. A Markov state model (MSM) of the ACBP folding reaction, constructed from over 30 ms of molecular dynamics trajectory data, predicts a complex network of metastable stables, residual unfolded-state structure, and kinetics consistent with experiment but no well-defined intermediate preceding the main folding barrier. Taken together, these experimental and simulation results suggest that the previously characterized fast kinetic phase is not due to formation of a barrier-limited intermediate but rather to a more heterogeneous and slow acquisition of unfolded-state structure.

Structure and dynamics of an unfolded protein examined by molecular dynamics simulation.

The accurate characterization of the structure and dynamics of proteins in disordered states is a difficult problem at the frontier of structural biology whose solution promises to further our understanding of protein folding and intrinsically disordered proteins. Molecular dynamics (MD) simulations have added considerably to our understanding of folded proteins, but the accuracy with which the force fields used in such simulations can describe disordered proteins is unclear. In this work, using a modern force field, we performed a 200 µs unrestrained MD simulation of the acid-unfolded state of an experimentally well-characterized protein, ACBP, to explore the extent to which state-of-the-art simulation can describe the structural and dynamical features of a disordered protein. By comparing the simulation results with the results of NMR experiments, we demonstrate that the simulation successfully captures important aspects of both the local and global structure. Our simulation was ~2 orders of magnitude longer than those in previous studies of unfolded proteins, a length sufficient to observe repeated formation and breaking of helical structure, which we found to occur on a multimicrosecond time scale. We observed one structural feature that formed but did not break during the simulation, highlighting the difficulty in sampling disordered states. Overall, however, our simulation results are in reasonable agreement with the experimental data, demonstrating that MD simulations can already be useful in describing disordered proteins. Finally, our direct calculation of certain NMR observables from the simulation provides new insight into the general relationship between structural features of disordered proteins and experimental NMR relaxation properties.

The interplay between transient α-helix formation and side chain rotamer distributions in disordered proteins probed by methyl chemical shifts.

The peptide backbones of disordered proteins are routinely characterized by NMR with respect to transient structure and dynamics. Little experimental information is, however, available about the side chain conformations and how structure in the backbone affects the side chains. Methyl chemical shifts can in principle report the conformations of aliphatic side chains in disordered proteins and in order to examine this two model systems were chosen: the acid denatured state of acyl-CoA binding protein (ACBP) and the intrinsically disordered activation domain of the activator for thyroid hormone and retinoid receptors (ACTR). We find that small differences in the methyl carbon chemical shifts due to the γ-gauche effect may provide information about the side chain rotamer distributions. However, the effects of neighboring residues on the methyl group chemical shifts obscure the direct observation of γ-gauche effect. To overcome this, we reference the chemical shifts to those in a more disordered state resulting in residue specific random coil chemical shifts. The (13)C secondary chemical shifts of the methyl groups of valine, leucine, and isoleucine show sequence specific effects, which allow a quantitative analysis of the ensemble of χ(2)-angles of especially leucine residues in disordered proteins. The changes in the rotamer distributions upon denaturation correlate to the changes upon helix induction by the co-solvent trifluoroethanol, suggesting that the side chain conformers are directly or indirectly related to formation of transient α-helices.